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We demonstrate that weak parametric interaction of a fundamental beam with its third harmonic 
field in Kerr media gives rise to a rich variety of families of non-fundamental (multi-humped) solitary 
waves. Making a comprehensive comparison between bifurcation phenomena for these families in bulk 
media and planar waveguides, we discover two novel types of soliton bifurcations and other interesting 
findings. The later includes (i) multi-humped solitary waves without even or odd symmetry and (ii) 
multi-humped solitary waves with large separation between their humps which, however, may not be 
viewed as bound states of several distinct one-humped solitons. 



I. INTRODUCTION AND MODEL 



Recently parametric wave mixing in Kerr media has at- 
tracted significant attention (see, e.g., Refs. Jl]-|j| where 
continuous wave (CW) interaction and parametric self- 
trapping were investigated). This theoretical activity has 
been backed up by experimental advances, e.g., a novel 
scheme for quasi-phase matched third harmonic genera- 
tion (THG) has been suggested However, in previ- 
ous works devoted to spatial solitary waves due to THG 
in planar waveguides, only families of fundamental self- 
trapped beams were considered. In the Letter ||, for 
example, where solitons due to third harmonic genera- 
tion were considered for a bulk medium geometry, higher- 
order modes were not discussed in detail. By higher-order 
we refer to beam shapes whose transverse intensity typi- 
cally has a multi-peaked structure and has higher energy 
than the single-peaked fundamental state. 

In this work we analyze in some detail the structure 
and bifurcation phenomena of higher-order bright spa- 
tially localized modes or 'solitons', which we do not use 
in a strict mathematical sense, since the models in ques- 
tion are not integrable. The spatial configuration is as- 
sumed to be such that there is a well defined propagation 
direction and the beams are localized in n transverse di- 
rections, with n — 1 representing a planar waveguide 
and n — 2 a, bulk medium. Specifically we study models 
representing (l+l)-dimcnsional and (2+l)-dimcnsional, 
weakly- anisotropic media with cubic nonlinearity, under 
the phase-matched condition that the fundamental wave 
is resonantly coupled to its third harmonic. This is a 
particular degenerate case of solitons supported by the 
four- wave mixing processes which is not completely 
understood yet in full generality. We assume that the in- 
teraction between the fundamental and third-harmonic 



waves includes the effects of parametric four-wave mix- 
ing, self-phase modulation, and cross-phase modulation. 

We closely follow the derivation procedure of Ref. [Q , 
assuming that the fundamental and the third-harmonic 
beams have the same linear polarization. The result is 
the following normalized (dimensionless) system of cou- 
pled equations: 
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where u and w are the fundamental and third harmon- 
ics, respectively. Also for the case of spatial beams 
V 2 = d 2 /dx 2 + d 2 /dy 2 in the (2+l)-dimensional case, 
or V 2 = d 2 /dx 2 in the (l-fl)D case. The parameter a 
measures the shift in the propagation constant, which is 
induced by the nonlinearity and is also dependent on the 
quality of wave- vector matching between the harmonics, 
with a = 3cr corresponding to exact matching, and z is 
the propagation distance. For the spatial soliton case 
the dimensionless parameter a is the ratio of the wave 
numbers of the harmonics and is equal to 3. Note that 
the system ([j]) may also describe temporal pulse prop- 
agation of resonantly interacting fundamental and third 
harmonics in optical fibers. For this physical situation 
V 2 = d 2 /dt 2 (t is the retarded time variable) and a is 
the absolute value of the ratio of second-order group ve- 
locity dispersions for the first and the third harmonics 
and may be any positive number. 

Radially symmetric stationary beams are described by 
real solutions u(r) and w(r) which are defined by the 
following system: 
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Here r = \/x 2 + y 2 and s = 1 for the (2+l)D case, 
whereas r = x and s = for the (1+1)D case. These 
localized solutions depend only on a single dimensionlcss 
parameter a. Analysis of the linear part of Eqs. (^) 
in the limit r — > oo shows that conventional bright soli- 
tons (with exponentially decaying tails) can exist only 
for a > 0. 

By 'bright symmetric' in the remainder of this paper 
we shall mean (2+l)D solutions when the intensity of 
each localized harmonic reaches a maximum at r = and 
(1+1)D solutions with u{r) — u{—r) and w(r) = w(— r). 
Thus, we shall only seek these solutions on the interval 
< r < oo even in the (1 + 1)D case. Note further 
that Eqs. (||) have odd symmetry, that is, if [u(r),w(r)] 
is a solution then so is [— u(r), —w(r)}. Thus all solutions 
must come in pairs, the second solution being simply a 
change in sign (a phase shift of n) of both harmonics. 
For the case s = 0, it is additionally possible to have 
solutions which are odd in both harmonics, or which are 
neither odd nor even. The latter type of solutions we 
shall refer to as being 'bright asymmetric'. In this paper 
we shall consider mainly the solitons of bright symmet- 
ric type, but shall also present some results about bright 
asymmetric (1+1)D solitons. Dark solitons (localized so- 
lutions with nonzero asymptotics) are out of the scope of 
this paper. 

The case a < also has physical meaning, but there 
one should expect to find quasisolitons, which are almost 
localized stationary states that have small periodic oscil- 
lations in their tails. See e.g. |||]U ^ or ^ nc definition, 
examples and for some issues surrounding them. Qua- 
sisolitons in this model will form the subject of another 
work. Here we shall concentrate almost exclusively on 
the case a > 0. 

Using a direct analogy with the theory of x^ 2 ' (second- 
harmonic generation or SHG) solitons (e.g. Ref. j9|), we 
start our analysis from the so-called cascading limit when 
a ^> 1. In this limit w ~ u 3 /(9a) and the equation 
for u approaches the cubic-quintic Nonlinear Schrodinger 
(NLS) equation. This scalar equation possesses a familiar 
class of fundamental bright solitons consisting of a sim- 
ple bell-shape [there are also higher-order families in the 
(2+l)D case]. These fundamental solitons can then be 
used as a starting point in the search for families of sta- 
tionary solutions using numerical methods. These meth- 
ods comprise a standard shooting method at fixed a, and 
a continuation method allied to solution using a relax- 
ation method for solving an appropriately defined two- 
point boundary- value problem for Eqs. (||). This latter 
technique can trace paths of solutions as a varies. We 
choose to characterize these solitons by the value of nor- 



malized total power which is one of the conserved quan- 
tities of the system ([!]) 



Pt, 



(\u\ 2 +3a\w\ 2 )dA. 



(3) 



Here the integration extends over the appropriate one or 
two-dimensional infinite cross-section A. 

The dependence of Ptot on a for a branch of solitons 
is usually, at least in the case of a fundamental solu- 
tion, closely related to its stability. A necessary con- 
dition for stability in the case of fundamental multi- 
component solitons is typically given by a generalized 
Vakhitov-Kolokolov (VK) criterion |^0| , which often also 
appears to be a sufficient condition for soliton stability 
(see, e.g., Refs. pT[|). However, the complexity of Eqs. (|l|) 
which, for example, possess collapse-type dynamics in 
the (2+l)D case, may lead to instability of fundamental 
solitons even for branches which are supposed to be sta- 
ble according to the VK criterion Thus, below we 
use power versus-a dependence only for classification of 
soliton families, leaving a full-scale stability analysis for 
future consideration. 



II. RESULTS FOR BULK MEDIA 



First we present the results for the (2+l)D case. Fig- 
ure [l] shows the variation of the normalized total power, 
-Ptot, with the normalized mismatch parameter a, for dif- 
ferent types of one- wave and two- wave localized solutions 
of the system ([[]) with s = 1 . The corresponding soliton 
profiles at various points along the presented branches 
are given in Figs. ||-^ 

The first class of localized solutions of the system (|l|) 
consists of one-frequency soliton families for the third 
harmonic Wq, which exist for all a > and represent 
scalar Kerr solitons described by the standard cubic NLS 
equation which follows from the second of Eqs. ([!]) at 
u = 0: 



'wo 



dr 2 



1 dwo 
r dr 



awo 



9^o = 0. 



(4) 



These single frequency solitons differ from each other by 
the number of zero crossings in their radial profiles so 
that we denote the corresponding families as Tq (no cross- 
ing), T\ (one crossing), Ti (two crossings), etc. Examples 
of one- wave solitons belonging to different Tj families are 
shown in Fig. |^. Note that the normalized power P to t is 
constant for each of the Tj families. For example, for the 
fundamental one-wave soliton family T (which are, in 
fact, Townes solitons of Ref. [jl2j) we have P to t ~ 11-70 
for all a > 0. 
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FIG. 1. (a) Bifurcation diagram for solitons (solid curves) 
and quasisolitons (dashed curves) of Eqs. (|l|) in the (2+l)D 
case. (b) Expanded portion of Fig. |l|(a) for the range 
< a < 10. Examples of solitons are shown in Figs. B-M. 
Bifurcation points of two- wave solitons from one- wave soliton 
families are shown by filled circles. The results related to qua- 
sisolitons are for stationary solutions with minimal amplitude 
of oscillatory tails; in that case, P is calculated for the soliton 
core only. 

The second class of solutions to Eqs. ([!]) are genuinely 
two-wave bright solitons, described by families of local- 
ized beams with coupled fundamental and third harmon- 
ics. The simplest way to obtain such solutions numeri- 
cally is to follow the two-wave soliton families from the 
cascading limit (large a) as a decreases. In this work 
we concentrate on the result of following the lowest or- 
der two-wave soliton branch whose profiles have a sim- 
ple one-hump form in the cascading limit. For this fam- 
ily, painstaking numerical continuation reveals a highly 
complex solution path involving restructuring of the soli- 
ton profile while the corresponding P(a) curve undergoes 
several loops (see Figs. |l| and ??). Inherent in each loop 
is a touch with one of the Ti families. Such a touch cor- 
responds to a transcritical bifurcation from the pure w- 
solution, and note [for example from Fig. @(c), (d) which 
correspond to points C and D on Fig. [j](aj] that the two 
different bifurcating branches have opposite signs of their 
it-component. The fact that these bifurcations take place 
further illustrates the severity of the restructuring of the 
soliton profiles that must take place; in the cascading 
limit the branch is approximately of pure u-type, whereas 



at each bifurcation with Tj it is composed of purely a 
third harmonic component w. Figures illustrate the 
complete restructuring process by depicting the soliton 
profiles in the vicinity of each bifurcation and turning 
point of the P(a) curve. Note finally that the two-wave 
soliton family also includes the simplest so-called self- 
similar [for which u(r) oc w (r)] solution (Fig. g, point M) 
at a = 1, see Ref. |iu| for the details and also Ref. [||, 
where its (l+l)-dimensional counterpart was also been 
considered. 



The position of the bifurcation point from the To 
branch can be approximately calculated analytically. 
Linearization of Eqs. (|l|) around the solution wo(r) gives 
the eigenvalue equation 



d 2 ui 
dr 2 



1 du\ 
r dr 



+ 2w 2 {r)ui = Ami, 



(5) 



together with appropriate boundary conditions. Bifur- 
cations occur when A = 1. This may also be viewed as 
the problem of existence of localized states in the po- 
tential U(r) = —2w 2 .(r) with eigenvalue A. Due to the 
lack of a closed form analytical expression for Wq(t), so- 
lutions of (^) may be approximated by feeding in the 
numerical data for wq or by analytical techniques based 
on a variational approximation. Using the latter, based 
on a simple exponential trial function, gives the result 
wo = V8a/3e~ r ^ . Substituting this into Eq. (||) and 
assuming a similar form of trial function for u\, one can 



use a Rayleigh-Ritz method to obtain a 



(var) 
bif 



105.8 



This agrees to within 2% with the numerical result 
a lif = 104.2. Calculation of bifurcation points along 
the higher-order Tj branches may in principle be carried 
out by the same method. However, this is less straight- 
forward technically because it requires the use of com- 
plicated forms of trial functions, and it is perhaps more 
instructive to rely on numerical detection of the bifurca- 
tion points. Symmetry arguments dictate that at each 
bifurcation point olu / there will exist two different bifur- 
cating solutions of Eq. (||): u\(r) and — u\(r). Moreover, 
each bifurcation is a transcritical and gives rise to a pair 
of two-wave branches (wq{t), ±£Mi(r)), where e is pro- 
portional to | a — otbif\. 

The third class of solutions to Eqs. (|l]) are the afore- 
mentioned quasisolitons which exist in the region of neg- 
ative a. We do not discuss quasi-solitons here in any 
detail. A full analysis will appear elsewhere. We simply 
make the comment that the branch SOP bifurcating from 
T3 can be continued up to the boundary a — separat- 
ing regular from quasi-solitons. On the other side of the 
boundary a similar quasi-soliton state can be found with 
tiny oscillations in its tail [see Fig. [j] and Fig. fj](t,s)]. 
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FIG. 2. Examples of different families of one-wave solitons. 
In all diagrams, the thick line corresponds to the 3 har- 
monic. 



FIG. 5. Examples of (2+l)D two-wave solitons. Labeling 
of examples corresponds to labeling of open circles in Fig. fl\ 





III. RESULTS FOR PLANAR WAVEGUIDES 



It is interesting to compare the (2+l)D results dis- 
cussed above with those for the corresponding (1+1)D 
case. The bifurcation diagram related to the (1+1)D case 
is presented in Fig. ^ and the corresponding examples of 
soliton profiles are given in Fig. ||-14. We now highlight 
how, together with many obvious differences in compari- 
son to the diagram for the (2+l)D case in Fig. [|, there are 
also some striking similarities as well. Note that in some 
respects the model for the (1+1)D case is simpler since 
the corresponding stationary system @ with s = does 
not depend explicitly on r and hence represents an au- 
tonomous dynamical system in four dimensions. Finding 
solitons is then reduced to finding homoclinic trajectories 
in this 4D phase space. 
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FIG. 8. (a) Bifurcation diagram for symmetric solitons 
(solid curves) and quasisolitons (dashed curve) of the (1+1)D 
version of Eqs. (|l|. (b) Expanded portion of Fig. |](a) in the 
range < a < 10, 40 < P < 65. Dotted curves emerging at 
zero correspond to integer multiples of the primary one-wave 
solution Si. Formally they represent multi-soliton states con- 
sisting of a concatenation of infinitely separated single soli- 
tons. Points at which branches of two- wave solitons terminate 
by 'bifurcating' from one of these multi-solitons are depicted 
by filled circles and all occur for a — 9. The inset to (a) and 
the jump (T — > N) depicted in (b) are explained in the text. 



The first class of (1+1)D localized waves of system (|l|) 
consists of one-frequency soliton families for the third 
harmonic wq, which exist for all a > and repre- 
sent scalar Kerr solitons described by the standard cubic 
(1+1)D NLS equation which follows from the second of 
Eqs. dj) at u = 0: 



d 2 w Q 
dx 2 



awo 



9w$ = 0. 



(6) 



It can be readily solved exactly giving the well-known 
unique single soliton solution: 



(7) 



w o(x) = ^-^sech(^/ax) , P to t = ^Va- 



in contrast to the (2+l)D case strictly speaking there 
are no other one-wave localized solutions. However, it 
will be helpful in what follows to consider formal multi- 
soliton states consisting of a different number of infinitely 
separated single solitons (0), families of which we de- 
note by S\ (single soliton), S2 (two solitons), S3 (three 
solitons), etc. In this work we are mainly interested 
in families with an odd number of separated solitons: 
SWfij i — 1)2,3,..., but we also investigate 'bifurca- 
tions' from S2. Note that, for i > 1, Si in fact denotes 
more than a single one-wave family, because each sin- 
gle pulse that is glued together can be either positive or 
negative. 

The second class of (l+l)-dimensional localized solu- 
tions of Eqs. (Q) consists of two-wave bright symmetric 
solitons and is described by families of localized beams 
with coupled fundamental and third harmonics. The sim- 
plest way to obtain the lowest order two-wave soliton 
family is again to continue numerically from solitons of 
the cascading limit (a » 1) given approximately by the 
expression: 



u(x) 



6 



y/l + B cosh 2x 



w pa u 3 /{9a), 



(8) 



where B = y/l + 16/a. The first expression for u(x) in 
Eq. (||) is the solution of the corresponding cubic-quintic 
NLS-type equation. 

The results of our numerical continuation from this 
limiting solution, upon decreasing a is that, like in the 
(2+l)D case, this branch also traces a convoluted path in 
the (P, a)-plane, involving four 'bifurcations' from one- 
wave soliton families (from the families S±, S3, S5, and 
57). As in the (2+l)D case, this branch connects to a self- 
similar solution at a = 1 (the point O in Fig. |8|(b)). In 
this case, the self-similar solution is expressible in closed 
analytical form as 



u(x) = a sech x, w(x) — b u(x), 



(9) 



where the parameter b is the real root of the cubic equa- 
tion 63fe 3 -3fe 2 + 176+1 = 0, and a 2 = 18/(186 2 + 36+ 1). 
However, it is here that the similarity with the (2 + 1)- 
case ends, as we shall now explain. 
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FIG. 9. Examples of (1+1)D two-wave and one-wave soli- 
tons. Labeling of all examples corresponds to the labeling of 
the open circles in Fig. pi 
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FIG. 12. Examples of (1+1)D two-wave solitons. Labeling 
is as for Fig. H. 
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FIG. 10. Examples of (1+1)D two-wave solitons. Labeling 
is as for Fig. ^. 
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FIG. 13. Examples of (1+1)D two-wave solitons. Labeling 
is as for Fig. ^. 
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FIG. 11. Examples of (1+1)D two-wave solitons. Labeling 
is as for Fig. ^. 
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FIG. 14. Examples of (1+1)D two-wave solitons, which are 
not directly linked to the two-wave solitons of the cascading 
limit. Labeling is as for Fig. ^. 



6 



First, let us try to motivate what is happening at each 
of the 'bifurcations' from Sj\ for which at first sight it 
seems remarkable that each one occurs precisely at a = 9. 
Standard bifurcation analysis (e.g. as in Ref. allows 
us to find the position of the single bifurcation point from 
the one-wave soliton family Si (Q) at a = 9.0 (point C 
in Fig. ||(a)). As in the (2+l)D case the bifurcation is a 
transcritical with one branch emerging to the left of the 
bifurcation point and one to the right. This structure is 
confirmed by the inset to Fig. ||(a) which shows that the 
branch emerging to the left undergoes a fold (at point 
B), so that on a larger scale both branches appear to 
bifurcate to the right. 

Now it seems that this 'local' bifurcation from Si 
causes a topological change in the four-dimensional phase 
space so that a global event must also happen at this 
parameter value. This global event is the possibility of 
gluing together several copies of the Si back to back 
and forming a new branch of solitons with several large 
peaks that bifurcate from a = 9. Phenomenologically 
this is similar to what happens in the SHG case when 
the parameter equivalent to a passes through 1 p5| , ^6[ . 
A key observation here is that in order to get a symmet- 
ric (even) solution, only an odd number of copies of the 
Si may be taken to form solitons in this way. As a con- 
venient short-hand for this global bifurcation of multi- 
peaked solutions at a — 9, we have refereed to it as a 
local 'bifurcation' from S2i+i where i — 1,2,3..., al- 
though this is strictly a misnomer. 

Numerical continuation beyond point G of Fig. ^(a) 
shows that two- wave soliton branch approaches a = 9.0 
from the left, where it bifurcates from the S3 asymptotic 
one- wave family that has alternative phase between each 
single-soliton component. However, we find that this is 
only one of a total of four symmetric two-wave solitons 
that come out of S3. There are 8 in total if you include 
the change of sign of both u and w. The second bifurcates 
to the left from the same (alternating phase) S3 family 
and differs only in that the first harmonic has the oppo- 
site sign. A representative of this branch, corresponding 
to point H in Fig. ||(a) is shown in Fig. [l0|(/i). The two 
other branches exist for a > 9 and bifurcate from the S3 
family where all peaks are in phase (positive) , and repre- 
sentatives are shown in Fig. |lj(u,w). With the increase 
of a (cascading limit) these complex multi-humped soli- 
tons keep their general structure intact, but become more 
localized. These two branches are not shown in the bi- 
furcation diagram (Fig. ||) but their P(a) curves lie very 
close to each other and to the T3 curve to the right of the 
bifurcation point. 

A similar bifurcation picture is observed at a = 9.0 for 
bifurcations from S5 and S7 one-wave families. However, 
because of the increase of the number in possible one- 
wave multi-soliton families themselves, the number of the 
corresponding bifurcated two-component branches also 
increases. For the even solitons considered in this work 
we have the following formula to calculate the number of 
two- wave sub-families bifurcating from one-wave Si fam- 



ily: Ni — 2^ +1 )/ 2 (double that if we count the change 
signs of u and v). For example, there are 16 branches 
that bifurcate from S7 branches which have P = 84 at 
a = 9. Note that in the bifurcation diagram of Fig. ||, 
in order to clutter, only branches directly linked to the 
cascading limit two-wave family are shown. Close to bi- 
furcation points, the third harmonic components of the 
depicted branches have neighboring humps of alternat- 
ing sign and first harmonic components have all humps 
of the same sign. Note that these branches all bifurcate 
to the left of a = 9. For the branches which bifurcate 
to the right not all third harmonic neighboring humps 
alternate in sign. 

It is important to note that none of the multi-hump 
soliton branches bifurcating to the left of a — 9 can be 
viewed as bound states of single partial solitons. Indeed, 
single one-hump solitons of Eqs. (Q) always have u and 
w components in-phase (of the same sign) for a < 9.0, 
whereas some of the individual humps of the of multi- 
hump structures bifurcating to the left from Si (i > 1) 
families have u and w components of different signs. To 
illustrate this point we show in Fig. 15 an enlarged bifur- 
cation diagram in the vicinity of a = 9 covering the first 
three families, Si, i = 1, 2, 3. Some of the corresponding 
examples of soliton profiles plotted at a — 8.6 are given 
in Fig. |l^. As they approach a = 9.0 the separation be- 
tween each individual hump (a 'partial soliton') increases 
and the state begins to approach a concatenation of single 
solitons with slightly overlapping tails. However, some of 
these partial solitons have out-of-phase u and w compo- 
nents and hence cannot exist on their own (i.e. without 
being in superposition with other 'partial' solitons). 
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FIG. 15. Bifurcation diagram from the first three 
one-component families Si,i = 1,2,3. Asymmetric family 
S2 is shown by a thick line. 
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FIG. 16. Examples ol the two-wave solitons close to bifur- 
cation point at q = 9. Weak component u(x) is enlarged in 
two bottom plots. Labeling of the profiles is in agreement 
with Fig. 0. 

Fig. ^ shows something even more striking - that 
there is also a 'bifurcation' from the 52 family. How- 
ever, the solitary waves that bifurcate from there are not 
bright symmetric but in fact are asymmetric solitons, see 
Fig. [l7|. Also at least one of these asymmetric solutions 
is born in a symmetry-breaking (pitchfork) bifurcation 
from one of the symmetric soliton branches (at the point 
O as , see Fig. |l^) . Thus there is a branch of asymmetric 
solitons which connects symmetric solitons with a branch 
of asymptotic antisymetric solitons (the S2 family) . We 
conjecture that there are similar asymmetric solitons that 
'bifurcate' from Sj at a = 9 for all even j. 
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FIG. 17. Examples of asymmetric solutions bifurcated 
from the family 52- Labeling of the profiles is in agreement 
with Fig. 

In contrast to the (2+l)D case, we have found no ex- 
amples (at least considering all bifurcations from S^i+i 
with 2i+l < 7) of two- wave solitons that survive down to 
a = where they might form a connection with branches 



of quasi-solitons existing for a < 0. Instead a repre- 
sentative branch coming from T7 bends abruptly (at R) 
at which point a increases through the point S until it 
reaches T at a rs 3.65, where another nonlocal bifurca- 
tion occurs. In this process, the third harmonic grad- 
ually forms a core with weakly separated wings. At T 
,The latter become completely separated one-wave soli- 
tons [see Fig. [l3](s,t)]. The solution at the point T can 
thus be viewed as a direct sum of two well-separated one- 
wave solitons and the soliton at point N. Beyond T we 
were unable to find any similar solutions. This non-trivial 
"jump" bifurcation is indicated by the vertical arrow in 
Fig. I 



IV. CONCLUSION 



In conclusion, we have investigated and classified 
higher-order soliton families and bifurcation phenomena 
due to resonant parametric interaction of a fundamental 
frequency wave with its third harmonic. 

In the case of (2+l)D solitons the picture is consis- 
tent with standard theories, albeit the branch we followed 
from the cascading limit connects several distinct soliton 
types in a non-trivial way. Also the structure of the sets 
of branches we found to approach the limit a = could 
do with further investigation, perhaps using singular per- 
turbation theory. The relation of these states for positive 
a to quasi-solitons for negative a will be addressed else- 
where. 

In contrast, in the (1+1)D case the bifurcation diagram 
is less clear-cut and we have found at least two novel fea- 
tures (i) the non-local bifurcation of multi-humped two- 
frequency solutions which are a consequence of the lo- 
cal bifurcation from the one-humped one-frequency soli- 
ton at a = 9, and (ii) the so-called jump bifurcation at 
the point T. The first of these is particularly intrigu- 
ing since not only are symmetric multi-humped states 
formed in this way, but also asymmetric ones. Also 
some of the multi-humped states cannot be viewed as 
bound states of several distinct one-humped states. The 
second novel bifurcation, the jump, appears related to, 
but not the same as, the so-called orbit-flip bifurcation 
p7[ . A dynamical-systems-theory explanation of these 
new bifurcation events, perhaps using the Lin-Sandstede 
method as in Ref. [Q, would be most interesting. 

Stability of the newly discovered soliton families re- 
mains an open question, especially for the (1+1)D case. 
Although usually higher-order soliton families are sub- 
ject to one of several types of instability, some exceptions 
are known (see, e.g., |Q) and thus a careful stability 
analysis is worth doing. The promise of detecting stable 
multi-hump solitons is real indeed because at least some 
of them cannot be viewed as bound states of two or more 
single (one-hump) solitons. For such bound state solitons 
of NLS-type system of equations, there is practically no 
hope of stability as shown e.g. in Ref. |19||. 
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